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Abstract. Elliptical galaxies are often assumed to be the 
primary source of heavy elements in the intracluster me- 
dium (ICM) , with the contribution of other morphological 
types being negligible. In this paper we argue that a role of 
spiral galaxies in the chemical evolution of the ICM is also 
important. This statement rests upon our recent calcula- 
tions of the heavy element loss from a disk galaxy through 
the hot steady-state galactic wind and dust grains expul- 
sion by stellar radiation pressure. This model reproduces 
main properties of our Galaxy and, being applied to galax- 
ies of various masses, explains the observed correlation be- 
tween spiral galaxy mass (luminosity) and metallicity. In 
our model this correlation develops as a result of the mass 
dependence of both loss mechanisms, in the sense that less 
massive galaxies lose metals more efficiently. We show that 
a typical disk galaxy is nearly as effective in enriching the 
ICM as an elliptical galaxy of the same mass. 

Having estimated the oxygen and iron loss from a 
single galaxy, we integrate them over the galactic mass 
spectrum. We show that the 'effective' loss (per unit lu- 
minosity) from spiral galaxies is comparable to the loss 
from ellipticals. The dominant role of early-type galaxies 
in rich clusters is caused by that they outnumber spirals. 
We present some arguments to this point, based on recent 
determinations of the ICM abundances, emphasizing the 
fact that the ratio of total iron mass to cluster luminosity 
does not depend on the fraction of cluster spirals in a wide 
range of the latter, contrary to what one might expect if 
spirals do not contribute into the ICM Z-abundance. 

Key v^rords: galaxies: clusters: general - galaxies: evolu- 
tion - galaxies: intergalactic medium - galaxies: spiral 



1. Introduction 

Clusters of galaxies are believed to be the only closed 
chemically evolving systems in the Universe, in the sense 
that they do not lose matter into the intercluster space. 
Inside them, however, violent processes of matter and en- 
ergy exchange between the galaxies and the intracluster 
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medium should occur. This conclusion comes from the 
fact that heavy elements, which are nearly equally dis- 
tributed between galaxies and the ICM, can only be pro- 
duced within galaxi es in t he course of 'normal' stellar evo- 
lution (e.g. Renzini 1997 ). 

Since in recent years the number of the ICM abun- 
dance determinations and their accuracy have increased 
significantly, the interest to the field of the chemical evo- 
lution of galaxy clusters has raised, too. Several mecha- 
nisms have been suggested to explain the metal transport 
from the galaxies to the ICM. Among them are enriched 
gas ejection during the merging of protogalactic fragments 



(protogalaxies) in the early Universe (Gnedin 199S), ram- 
pressure stripping of enriched gas both from spiral and el- 
liptical galaxies (Gunn & Gott 1972; Himmes & Biermann 
1980; Fukumoto & Ikeuchi 1996), and the hot galactic 



wind from ellipticals (Matteucci & Vettolani 1988; David 
et al. |1991| ; Matteucci & Gibson |l99|; Renzini |1997| ; Gib- 
son & Matteucci 1997[ and others). 

From these, the last mechanism draws more and more 
attention as an ultimate explanation of the ICM Z-abun- 
dance. Moreover, as Arnaud et al. (1992) pointed out, it 
was the hypothesis about the galactic wind from ellipti- 
cals, which led to prediction of heavy element existence in 
the ICM (Larson |1974| ; Larson & Dinerstein |l975[ ). This 
hypothesis arises from the assumption that ellipticals be- 
gan their history with the strong starburst followed by 
multiple supernova events. These explosions after a very 
short time expelled all the interstellar matter into the 
ICM, having provided it with heavy elements and thermal 
energy. At the same time, hot galactic wind is thought to 
be responsible for establishing the correlation of luminos- 
ity (mass) and metallicity, observed in ellipticals. 

As several theoretical studies of this process have 
shown, mass loss from ellipticals is effective enough to ex- 
plain Z-abundance in the ICM, though this mechanism, 
according to some of these studies, requires the IMF of 
high z galaxies to be richer in massive stars than today 
(Matteuc ci fc V ettolani [1988| ; David et al. |1991| Matteucci 
& Gib son [1995 ; Okazaki et al. 1993). Futhermore, Arnaud 
et al. ( 1992| ) found observational confirmation of the ellip- 
ticals being the dominant source of intracluster metals in 
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the apparent increase of the iron mass in galaxy clusters 
with the growth of the combined E+SO luminosity. From 
this they concluded that the role of spiral galaxies in the 
ICM enrichment is negligible. 

This possibly explains why spiral galaxies lack at- 
tention as a potential source of the intracluster metals, 
though, as Fukumoto & Ikeuchi (1996) argued, a single 
spiral galaxy could be as effective in the ICM enrich- 
ment as an elliptical one, and correlation of iron mass 
with the combined luminosity of ellipticals and lenticulars 
is caused by that early- type galaxies outnumber spirals. 
Furthermore, recently Kauffmann & Chariot ( 199§| ) sug- 
gested that metals in the ICM had been ejected from the 
primordial population of spirals that have later merged to 
form present-day early-type galaxies. In their model, spi- 
ral galaxies are not the dominant, but the only source of 
heavy elements in hot intracluster gas. 

Recently we also argued that disk galaxies suffer sig- 
nificant Z-loss (Shustov et al. 1997; hereafter Paper I). 



Tsujimoto et al. (1995; SN la). Each la Supernova is as- 
sumed to produce O.QMq of iron. Time delay for these SNe 
is taken to be equal to 3 x 10^ years (Tutukov & Yungel- 
son 1994). We use Salpeter-like initial mass function with 



We considered two processes leading to heavy element ex- 
pulsion from disk galaxies: hot galactic wind and dust ex- 
pelling by the stellar radiation pressure. It was shown that 
these processes are rather effective and may result in the 
luminosity-metallicity correlation observed in disk galax- 
ies. In this paper we extend these calculations to the in- 
tegration of heavy-element production in disk galaxies of 
various masses over the galactic mass function, thus as- 
sessing Z-loss from all disk galaxies in a typical cluster. 

In Section 2 we briefly describe the chemical evolution 
model used for computation of heavy element loss from a 
disk galaxy. In Section 3 results of integration of oxygen 
and iron losses over the galactic mass function are pre- 
sented and compared with the analogous data for ellipti- 
cal galaxies. They are discussed in Section 4. A conclusion 
follows that a disk galaxy enriches the ICM as effective as 
an elliptical one of the same mass. 

2. Chemical evolution model 

Model of galactic evolution used here is outlined in Paper I 
and Wiebe et al. ( 1998). Its primary feature is that in addi- 



tion to the conventional mass conservation equations for a 



total gas mass and for a particular element (Tinsley 1980) 
we solve the energy equation, based on a balance between 
the turbulent energy input from supernova explosions and 
its dissipation in cloud-cloud collisions. In Paper I we have 
successfully used this model to show that disk galaxies suf- 
f er sig nificant loss of heavy elements, while in Wiebe at al. 
(1998) we showed that this model is good in reproducing 
main characteristics of our Galaxy. 

As in Wiebe et al. ( 1998| ), here we model the evolution 
of oxygen and iron. For the oxygen production wc interpo- 
late metallicity-dependent yields of Maeder ( |1992D . Iron is 
assumed to originate both in core-collapse and la super- 
novae. Its production is taken into account according to 
Thielemann, Nomoto & Hashimoto (|l995| ; SN Ih+U) and 



limits 0.1 and IOOMq. 

Within the context of the present paper it is impor- 
tant to remind how we take into account heavy element 
loss from disk galaxies. We consider two processes eject- 
ing Z into the ICM. They are hot galactic wind and dust 
expelling by the stellar radiation pressure. Heavy element 
loss rates due to both processes depend a on galaxy mass, 
in the sense that low-mass galaxies lose their metals more 
effectively. 

Hot galactic wind forms from merged galactic foun- 
tains produced by multiple supernova explosions in OB- 
assoc iations (Mac Low et al. 1989 ; Nor man fc Ikeuchi 
1989t Igumentschev, Shustov, & Tutukov [19901) . Only a 
fraction /esc of this wind is lost into the circumgalactic 
space. Combining observational and theoretical results, we 
derived in Paper I linear dependence of /esc on the galaxy 
mass. This dependence is based on two theoretical and 
several observational points. According to Igumentschev 
et al. ( 199C| ), /esc ^ 1 for galaxies with masses greater 
than 10^"^ Mq. De Young & Gallagher (|1990D showed that 



in a galaxy of 1.4 x 1O^M0 /esc < 0.6. In Paper I we 
supplied these theoretical limits with several intermediate 
observational points, assuming that the surface filling fac- 
tor of supershells produced by the multiple explosions can 
be a natural estimate for /esc- Resultant dependence of 
/esc on the logarithm of the galaxy mass is well fitted by 
a straight line (see Paper I for details). 

The second process leading to Z-\oss from a disk 
galaxy is dust expelling by the stellar radiation pressure. 
Existence of extended dust halo around disk gala xies is 
demonstrated in a nu mber of observations (Zaritski 1994; 
Howk & Savage 1997 ). Using numerical model elaborated 
by Shustov & Wiebe ( 1995 ), we estimated dust loss de- 
pendence on a galaxy mass. It turns out that effectiveness 
of dust expelling decreases with galactic mass, again ap- 
proaching nearly zero for galaxies with masses ~ lO^^M©. 

A new important feature is included in the present 
version of our model. In its previous version, described in 
Wiebe et al. ( 199^ ), the galaxy forms at once that leads to 
high star formation rate in the early galaxy life. To make 
our model more realistic, in its present version we change 
the initial conditions by adding the initial accretion phase. 

In a collapsing (protogalactic) cloud, the rate of matter 
accretion onto the cloud core depends only on the cloud 
initial temperature (e.g. Massevitch & Tutukov 1988| ). We 
assume in our standard model that the accretion rate is 
lOOM© yr^^. That corresponds to the accretion time 2 x 
10^ years. During this time mass of the Galaxy grows with 
constant rate; when it equals to Mq, the accretion stops. 
This makes the starburst lower and smoother. 
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Fig. 1. Redshift (time) dependence of oxygen (a) and iron (b) masses ejected from a galaxy with Mq = 5 • 10^^M( 
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Fig. 2. Dependence of theoretical oxygen (a) and iron (b) ejected mass on the mass of a galaxy for spirals (this paper) 
and ellipticals (previous investigations). 



3. Results 

We perform calculations as follows. Using our Galaxy 
as a standard, we find proper values of free parameters 
to reproduce its main properties, including dependence 
of metallicity on the height above the galactic disk and 
[0/Fe]-[Fe/H] correlation. After that, adopted values of 
parameters are kept fixed and used further to model galax- 
ies of different masses. 

We investigate the evolution of oxygen and iron ejec- 
tion from galaxies with masses between 10^° and 5 • 
IQ-'^^Mq. Less massive galaxies also eject heavy elements 
into the ICM, however, as we noted in Wiebe et al. (1998D, 



our model probably does not work for them, mainly be- 
cause the assumption of the self-regulated star formation 
is not valid in small systems. In order to compare our 
results with existing estimates of heavy element ejection 
from elliptical galaxies we compute two models: Salpeter 
IMF {x = 1.35) and flat IMF (x = 0.95). 



Shown in Fig. n^is a time dependence of ejected oxygen 
and iron mass from the galaxy with Mq = 5 • IO^^Mq for 
the two different cases (throughout the paper the standard 
cosmological model with h = 0.5 is adopted). Introducing 
in the model the flat IMF with higher proportion of mas- 
sive star increases both oxygen and iron production nearly 
by an order of magnitude. 

In Fig. 0a our results for oxygen loss (integrated over 
time to i = 13-10^ years) are presented in comparison with 



predictions of David et al. (1991) and Matteucci & Gibson 
(1995). Fig. &) depicts results for iron ejection according 
to the same papers. In all cases upper curve corresponds to 
the flat IMF, and lower curve corresponds to the 'ordinary' 
IMF (x = 1.35 for this paper and Matteucci & Gibson; 
X — 2.0 for David et al.). 

As it is demonstrated by Fig. 0, disk galaxies with 
masses less than 5 ■ IO^^Mq are as effective in enriching 
the ICM with O and Fe as ellipticals. Specific feature of 
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Fig. 3. Integrated oxygen (a) and iron (b) ejected mass vs. redshift for different choices of stehar and galactic mass 
function. 



our model is steep decrease in the heavy element ejec- 
tion rate for more massive galaxies. This is caused by the 
adopted form of dependence of the ejection efficiency on 
the galactic mass, as explained above. The similar decrease 
of the heavy element ejection rate for massive d isk ga lax- 
ies was argued also by Kauffmann & Chariot ( 1998 ). In 
their model, galaxies with circular velocities greater than 
250 km s""'^ (i.e., slightly more massive than the Galaxy) 
eject only 20 per cent of the intracluster metals. 

To assess the overall production of O and Fe in disk 
galaxies one has to integrate mass of the element ejected 
from a single galaxy over the galactic mass function. It 
is now widely assumed that the present galactic luminos- 
ity function (LF) can be described by Schechter (1976) 
law with reasonable accuracy. We take Schechter function 
in its original form, assuming that the spiral luminosity 
function is proportional to it (universality of LF for all 
morphological types was recently confirmed by Andreon 



1998 and Marzke et al. 199^) 


N = fn^ -— exp 





J-J^ 



(1) 



with limits li and I2 being the minimum and maximum 
luminosities, expressed in units of L*. Here / is a number 
fraction of spirals and n, is cluster richness, as defined by 
Schechter ( 1976| ). One of the quantities computed in our 
model is final galaxy luminosity. Comparing it with the 
galaxy mass, we find that in our model the present-day 
mass-luminosity relation fits the power law 

L = 1.3A-f"-^. 

Substituting it in (|l|), we obtain the present-day mass 
function 



N = mn.n^ 



0.9Q-0.1 



exp 



M 



0.9' 



-1^ 



where mi are mass limits mentioned above and Af* = 
(L*/1.3)^/°'^. Now we may write the formula for the mass 
of this element ejected by all galaxies in the cluster 



M^°* 



X exp 



0.9M 



M 



0.9Q-0.1 



M 



Mzd 



M 



(2) 



where Afz is the mass of heavy elements ejected by the 
galaxy of the mass M . Using this function we can com- 
pute the evolution of the intergalactic oxygen and iron 
produced in disk galaxies. 

Results of integration are presented in Fig. where 
oxygen and iron yields are shown as functions of redshift. 
They are computed with a = —1.3 and L* correspond- 
ing to absolute magnitude M-q — —22.5 (rich cluster case 
from Matteucci & Gibson 1995 ). As it was the case for a 
single galaxy, again the amount of the ejected element is 
increased by an order of magnitude for the flat IMF. Note 
that to convert numbers on y-axis into the total mass of 
the element ejected in the ICM by all disk galaxies in the 
cluster, one should multiply them by /n*. 

To compare our results with the predictions about O 
and Fe ejection from elliptical galaxies presented by Mat- 
teucci & Gibson ( 1995 ) we computed two models corre- 
sponding to their cases of a rich and poor cluster. Pa- 
rameters of these models and their results are listed in 
Table H. It is obvious that elliptical galaxies are the domi- 
nant source of heavy metals in rich clusters, while in poor 
clusters, rich in spirals, the latter ones catch up with el- 
lipticals in iron production and exceed them in oxygen 
production. We shall return to this point in the next sec- 
tion. 

To check dependence of our results on the adopted 
parameters of the LF we computed oxyge n and iron yields 
with a and L* from Marzke et al. ( 1998 ). Evolution of O 
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Fig. 4. Rate of oxygen (a) and iron (b) ejection for our 
Galaxy as a function of time. 

and Fe ejected mass for this LF is also shown in Fig. g|. In 
this case production of both iron and oxygen is decreased 
by an order of magnitude. However, predicted mass of O 
and Fe ejected into the ICM by ellipticals would decrease, 
too, so this cannot change results of comparison for the 
two morphological types. 

4. Discussion 

Hypothesis on a hot galactic wind associated with the 
burst of star formation in elliptical galaxies, being satisfac- 
tory in explaining temperature and abundance of the ICM, 
contradicts to an apparent lack of bright ellipticals in the 
young Universe, found in variou s deep field surveys (Kauff- 
mann. Chariot, & White |l996t Zepf |l997D . A number of 
effects has been proposed to explain this fact, intragalac- 
tic and intergalactic extinction being among them. An- 
othe r exp lanation was suggested by Kauffmann & Char- 
lot ( 1998 ). They argued that ellipticals did not form in 
a single burst of star formation but formed later during 
merging of disk galaxies. That is why we do not see their 
bright early stages. 

In this paper we do not intend to go into details of the 
elliptical galaxies formation. Our goal is to investigate the 
ICM-enriching role of a single non-interacting disk galaxy. 
However, our model has to deal with the problem of the 
initial starburst, too. Shown in Fig. H is the evolution of 
the rate of heavy element expulsion for our Galaxy (con- 
sidered as a standard). It is obvious that, though the rate 
is high enough even at the present epoch, the bulk of Z- 
loss occurs at early stages of galaxy's life as a result of 
intense star formation and related high supernova rate. 
Thus, this initial starburst is an important factor deter- 
mining the effectiveness of heavy element los s from disk 
galaxies. As we pointed out in Wiebe et al. ( 1998 ), it is 
impossible to avoid the initial starburst within the frame- 



work of our model without breaking consistency between 
the current theoretical parameters of disk galaxies and ob- 
servational data. Similar phase of the high star formation 
rate (though not so intense as in our model) is obtained, 
e.g., in hydrodynamical modelling of our Galaxy by Rai- 



teri et al. ( |1996D and Samland et al. ( |1997| ). 

There is a growing body of evidence that the star for- 
mation bursts did occur in most galaxies but is hidden 
from optica l surveys by the d ust ab sorption (see e.g. Blain 
et al. |1998|). In Wiebe et al. (|1998D we argued that at the 



time of the starburst a typical galaxy is able to accumu- 
late enough dust to screen out stellar radiation. In the 
most favorable conditions absorption can amount to sev- 
eral magnitudes. 

Another factor, strongly affecting our results, is the 
used method of accounting for the dependence of Z-loss 
efficiency on the galactic mass. The 2D hydrodynamical 
mode ls of a blow-out of the galactic disk (Igumentschev 
et al. 1990|), being a useful instrument for estimate, can- 
not be considered as a full treatment of gas dynamics of 
heavy element loss from the Galaxy. Note, however, that 
in dense gaseous disks, where all processes governing the 
self-regulated star formation (shock waves, ionizing radi- 
ation, etc.) have relatively short distance scales, the star 
formation efficiency is determined by the local gas param- 
eters and does not depend on global characteristics of a 
galaxy. This means that the final gas metallicity within 
the framework of our model does not depend on the mass 
of a galaxy. This fact is well illustrated by Table || where 
final metallicities are given with and without taking into 
account mass-dependent Z-loss. Thus, we have no way 
of reproducing the observed mass-metallicity correlation 
other than mass-dependent galactic wind and dust grain 
expelling. Keeping in mind that we are able to obtain the 
above correlation in our model (probably, underestimat- 
ing Z-loss from low-mass galaxies), we may say that this 
model is reasonably accurate. 

Apparently, we have a contradiction with the results of 
Arnaud et al. (1992), who argued for the hot wind from el- 
lipticals as the only significant mechanism of ICM enrich- 
ment with heavy elements. According to our data, both 
elliptical and disk galaxies give comparable contributions 
into the chemical evolution of the ICM. One of the basic 



premises of Arnaud et al. (1992) paper was the correlation 
on the intracluster gas mass - luminosity diagram for el- 
lipticals and the scatter on the same diagram for spirals. 
One may conclude from this fact that ellipticals, in gen- 
eral, play the dominant role in the intracluster medium 
evolution, while the spirals do not. However, as the theo- 
retical modelling of elliptical galaxies evolution show, all 
the gas ejected from the cluster ellipticals (in a typical 
cluster) cannot account for its amount in the ICM. So, 
the bulk of the intraclust er gas should have primordial 
origin (see Okazaki et al. 1993 for discussion), and the 
lack of correlation between the intracluster gas mass and 
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Table 1. Predicted mass of elements ejected into the ICM from all cluster E-fSO galaxies and spirals. 



Cluster 


a 


M5 


n. 


/ 


E+SO 
O Fe 


Spirals 
O Fe 


Rich 
Poor 


-1.3 
-1.3 


-22.5 
-22.0 


115 
20 


0.2 
0.7 


1.10(11) 4.04(10) 
3.77(9) 1.19(10) 


3.9(10) 7.9(9) 
1.9(10) 3.8(9) 
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Iron (a) and silicon (b) masses in the galaxy clusters vs. partial luminosity of galaxies of various morphological 



Table 2. Present-day oxygen abundance in open and 
closed models. 



Mass, 


Oxygen abundance [0/H] 


log(MG/M0) 


Open model 


Closed model 




Stars Gas 


Stars Gas 


10 


-0.4 -0.3 


0.2 0.1 


11 


-0.1 -0.1 


0.2 0.1 


12 


0.1 0.1 


0.2 0.1 



spiral luminosity, being of great interest on its own, has 
only a little bearing on the problem investigated. 

In 1992 Arnaud et al. had only six clusters with 
known metallicities and morphological populations. Re- 
cently, new data about the ICM chemi cal co mposition 
have become available (Fukazawa et al. 199S). Here we 
combine them with morphological data compiled by Ar- 
naud et al. ( [1992[ ) to understand if our results disagree 
with the current observations. 

In Fig. @ we show iron and silicon mass in the clusters 
versus partial luminosity of spirals, ellipticals, and lentlc- 
ulars (all data needed to compute these quantities were 
taken from Fukazawa et al. and Arnaud et al.). These 
two plots reveal the surprising fact: both iron and silicon 
masses follow nearly the same path for all morphological 
types, with the only exception of A426 (Perseus cluster). 
It is noteworthy that according to Andreon ( |1994 ) the 
number of spirals in this cluster is underestimated by a 
factor of 2 or 3. 



However, one might argue that the fact that iron (or 
silicon) mass in the cluster correlates with a partial lu- 
minosity of an y mor phological type is not so meaningful. 
Rcnzini et al. ( 1993 ) suggested that it is the ratio of iron 
mass to cluster luminosity (IMLR) that docs make sense 
in studying the ICM chemical evolution. If we plot fraction 
of spirals versus total luminosity for galaxies from Arnaud 
et al. ( 1992 ) sample (Fig. ^, we see that in the range of 
luminosities from 10^^ to lO^'^L© fraction of spirals drops 
from 0.5 to approximately 0.1. One might expect that the 
IMLR should decrease in spiral-rich clusters because of 
decreasing number of metal sources. 

To check if this is the case, we plot in Fig. R the ICM 
iron mass to light ratio (IMLR) and silicon mass to light 
ratio (SMLR) versus fraction of spirals in galaxy clusters. 
In both cases there is an apparent correlation between 
these two values. However, it rests mainly upon the cluster 
A426 which we have already mentioned. So, if we drop this 
cluster from consideration. Figs g| and reveal that, first, 
iron and silicon masses in the galaxy clusters increase with 
the luminosity of any morphological type. Second, there 
is no clear correlation between the IMLR (or SMLR) and 
a fraction of spirals in the cluster. 



5. Conclusion 

In this paper we investigate the role played by disk galax- 
ies in the chemical evolution of the ICM. We show that if 
a mass of the galaxy does not exceed ~ 5 • lO^^M© it con- 
tributes nearly the same amount of metals into the ICM as 



0.04 



0.03 



0.02 



0.01 




Wiebe et al.: Chemical enrichment of the intracluster medium 



0.04 



0.0 0.1 0.2 0.3 0.4 0.5 
Fraction of spirals 



0.03 



w 



i0.02 



0.01 



0.00 




0.0 0.1 0.2 0.3 0.4 0.5 
Fraction of spirals 



Fig. 7. Iron mass-to-light ratio (a) and silicon mass-to-light ratio (b) for a sample of galaxy clusters. 
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Fig. 6. Fraction of spirals as a function of total cluster 
luminosity. 

an elliptical galaxy of comparable mass. However, we have 
found that massive spirals are able to retain their heavy 
elements and thus do not participate in the ICM enrich- 
ment. Integration of the heavy element yields of galaxies 
over the galactic mass spectrum shows that while ellipti- 
cals do play the determining role in the chemical evolution 
of rich clusters, spirals can be significant source of Z in 
poor clusters. 

We should mention briefly some moments which are 
important for future improvement of our model: there are 
evidences that the galactic mass function changed in time 
(e.g. because of coalescense of early galaxies); new impor- 
tant data on abundance of other chemical elements are 
expected from space experiments; dynamical interaction 
of gas and galaxies in the clusters seems to be important 
and it should be investigated in more details. 
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